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We present a study of orbital angular momentum transfer from pump to down-converted beams 
in a type-II Optical Parametric Oscillator. Cavity and anisotropy effects are investigated and 
demostrated to play a central role in the transverse mode dynamics. While the idler beam can 
oscillate in a Laguerre-Gauss mode, the crystal birefringence induces an astigmatic effect in the 
signal beam that prevents the resonance of such mode. 

PACS numbers: 42.50.Dv; 42.65.Yj; 42.30.-d; 42.50.Lc 
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I. INTRODUCTION 

Early experiments have shown that circularly polar- 
ized light carries angular mometuni 0. In a quantum 
description of light, this angular momentum is associated 
with the spin of the photon. More recently a significant 
attention has been given to the study of the orbital angu- 
lar momentum of light, associated with phase singulari- 
ties in the wavefront. In a paraxial description of wave 
propagation, it is found that Laguerre-Gaussian beams 
carry orbital angular momentum. Such beams can be 
experimentally produced either through astigmatic mode 
conversion with cylindrical lenses 0, [4] or through holo- 
graphic techniques |j. |5j. |(| . While the holographic tech- 
niques are simpler, the astigmatic mode converters can 
operate with high power. 

It has been demonstrated that parametric down con- 
version with free propagating beams conserves the Or- 
bital Angular Momentum (OAM) of the pump beam. In 
the quantum optical domain this was observed through 
intensity correlation measurements, where twin-photon 
coincidences were obtained whenever the conservation 
condition was fulfilled 0. In the classical counterpart, a 
stimulating beam was introduced and parametric amplifi- 
cation was shown to be conditioned to OAM conservation 
[j| . Similar effects have already been observed in second 
harmonic generation 0, ITol ] . 

So far, little attention has been given to OAM con- 
servation in intracavity nonlinear coupling. Many stud- 
ies have been done with transverse multimode OPOs, 
showing interesting possibilities in pattern formation and 
quantum images for cavities with degenerate transverse 
modes, like planar 0,^3 an( i spherical cavities [pj.ll4|. 
Moreover, experiments have shown pattern formation in 
confocal and concentric ^(| cavities, and oscillation 
in modes with higher order than the fundamental are 
common in many different experiments [TtI Il8| . 

Apart from some theoretical studies on generation of 
phase singularities with nonlinear optical effects [T^. l2fj| , 
only a few experimental results have been published in 
this subject, and to our knowledge, there is no result 
showing the necessary conditions for intracavity OAM 



transfer from the pump to the down-converted beams. 
In the present work, we study the OAM transfer in a 
non-degenetrate, type-II Optical Parametric Oscillator 
(OPO). We show the conditions that must be satisfied 
for the OAM transfer, allowing one of the down-converted 
beams to oscillate with the same phase singularity of the 
pump beam. As we shall see, the astigmatism caused 
by the crystal birrefringence plays a central role in the 
selection of the beam oscillating in the Laguerre-Gauss 
mode. 



II. EXPERIMENTAL SETUP 

The experimental setup is shown in Figd The OPO 
is made by two spherical mirrors Ml and M2, with equal 
curvatures R m — 13mm. Inside the cavity, we have a 
KTP crystal (by Cristal Laser) 10 mm long, cut for non- 
critical phase matching in 532 nm to 1064 nm down con- 
version at room temperature. In this case, the crystallo- 
graphic axes (x, y, z) of the cristal are oriented as follows. 
The z axis of the crystal is vertically oriented while the 
propagation direction lies on the horizontal plane (xy). 
The x axis forms an angle cf> — 23.5° with respect to the 
propagation direction. 

The mirrors have high reflectance for the infrared (R 
= 99.8 % @ 1064 nm), and a small transmition at the 
pump wavelength (R = 92 % @ 532 nm). Crystal losses 
in the infrared comes mainly from surface reflection, re- 
duced by anti-reflective coating (R = 0.1 %), since crystal 
absorption at this wavelength is small (0.05 %). For the 
pump, we have reflection losses (R = 0.5 %) and crystal 
absorption, increased by gray-tracking effects [21]. The 
cavity length is controlled by a piezo actuator on the 
mirror, and the cavity is kept nearly confocal, in order to 
help the alignment and reduce the consequences of the 
walk-off. 

The OPO is pumped by the second harmonic of a 
Nd:YAG laser (Lightwave 142). This laser generates a 
TEMoo gaussian beam, that is converted to a nearly 
Hermite-Gauss TEM i beam . With a telescope 
formed by two cylindrical lenses, we implemented a mode 
converter that produces a Laguerre-Gauss beam 0, with 
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FIG. 1: Setup for the study of the phase singularities in the 
output of the type-II triply resonant OPO. The mode con- 
verter introduced in the pump beam is not shown 



a good cylindrical simetry for the intensity and a phase 
singularity in the center. This phase singularity was ev- 
idenced by the self interference pattern obtained in a 
Michelson interferometer. In Fig[21 we show the trans- 
verse profile and interference pattern of the beam used 
to pump the OPO. The resulting pump power is 60 mW. 
The beam is horizontally polarized, and mode matched 
to the cavity with the help of coated lenses. Although 




FIG. 2: a) Transverse profile of the pump; b) Interference 
pattern showing the topological defects characteristic of phase 
singularities 

the mirrors were high reflecting at 1064 nm, the output 
power coming out from the cavity through M2 can be 
detected by a PIN photodiode. The green light coming 
from the cavity is filtered by a dichroic mirror (DM), and 
detected by an amplified Si photodetector (D<y). The in- 
frared light is detected by a PIN InGaAs photodiode Dm 
(ETX-300, from Epitaxx), that samples part of the out- 
put beam that is reflected by a beam splitter (BS). 

In the output, signal and idler beams are separated by 
a polarizing cube (PBS). Adopting the usual convention 



in type-II OPOs, the idler beam polarization is horizon- 
tal, and the signal beam has a vertical polarization, al- 
ligned to the crystal z axis. Each down converted beam 
is sent into a Michelson interferometer made by a nonpo- 
larizing 50/50 beam splitter (BS) and two flat mirrors, 
in order to produce interference fringes that can reveal 
the existence of a phase singularity. The two outputs are 
recombined in another polarizign cube (PBS) and sent 
onto a CCD camera, that is used to register either the 
interference pattern or the intensity profile of the beam. 

The output power of the pump and infrared beams is 
measured as the cavity length is scanned. The corre- 
sponding resonance peaks are shown in FigO A wide 
peak is obtained for the pump, over which narrow dips 
appear due to the pump depletion in different oscillation 
regimes. The resonance peaks for the infrared are also 
shown in FigO They coincide with the depletion dips in 
the pump resonance. Expanding the curve, we can ob- 
serve that the depiction dips have a parabolic shape, in 
good agreement with the depletion expected for a triply 
resonant OPO 23]. From the finesse of the resonance 
peak for the pump, we measure 29 % of internal losses in 
the cavity. For signal and idler modes, the fitting of the 
parabolic depletion gives a value of 1 % for the infrared 
losses. The threshold power for parametric oscillation is 
around 20 mW. The OPO could be held in the oscillating 
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FIG. 3: Resonance peak for the pump beam, showing the 
oscillation of the OPO. The infrared peaks are labeled in order 
to identify the different images shown in FigHI 

peaks, and it would rest with a continuous output for as 
long as 10 minutes. In this situation, we registered the 
output image of signal and idler beams, as well as their 
self interference patterns. These images are shown in 
FiglU They are lebeled in correspondence to the oscilla- 
tion peaks shown in FigEI I n images 1 and 4, the output 
intensity in the idler is the one of a Laguerre-Gauss beam. 
The corresponding interference patterns show the topo- 
logical defects in the center of the Laguerre-Gauss beam 
characteristic of phase singularities. In this situation, the 
idler beam carries the orbital angular momentum of the 
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FIG. 4: a) Intensity patterns for signal (right) and idler (left) 
beams labeled in correspondence to the infrared peaks shown 
in FigE| b) Self interference patterns showing the presence 
or not of phase singularities. 



III. PARAXIAL WAVES IN A 
BIRREFRINGENT CRYSTAL 

In a type-II down conversion, we use the crystal bir- 
rcfringence to achieve the desired phase matching con- 
dition [2{J |2(J . It was shown by many authors that the 
paraxial equation, from which we can derive the propaga- 
tion modes of a beam in free space or isotropic medium, 
will change when we work with an anisotropic medium 
Here we will extend the study of Fleck and Feit 
|29j of paraxial propagation in uniaxial crystals to the 
biaxial case, adapting their description to the case of a 
crystal inside a cavity. Our aim is to reduce the wave 
equations to the paraxial wave equations that define the 
Hermite-Gauss modes coupled to the resonances of a lin- 
ear cavity. 

Let us define the crystalographic axis as (x, y, z). The 
KTP crystal used in our experiment is a quasi-uniaxial 
one with n x ~ n y ^ n z , where n x ( y ,z) is the refraction 
index for x(y, z) polarized light. Since the displacement 
vector D satisfies V • D = 0, we can write the wave 
equation for the electric field E, derived from Maxwell's 
equations, as 



V 2 E- V V • E 



V D 



k 2 £ -~E 



0, (1) 



down converted pump photons. In image 2, the shape 
of the idler beam is intermediary between a first order 
Laguerre-Gauss and a diagonal first order Hermite-Gauss 
modes. A vortex can still be observed through the inter- 
ference fringes. In both cases, the signal beam remains in 
the fundamental gaussian mode. Following the Poincare- 
sphere representation proposed in Ref . [24[ , we can look at 
the idler mode shown in image 2 as an orbital equivalent 
of an elliptical polarization. 

An interesting effect appears in image 3. In this situ- 
ation, the signal beam oscillates in the transverse mode 
with higher order, but with no angular momentum. The 
output is a pure Hermite-Gauss TEMoi mode, vertically 
oriented, while the idler remains in the fundamental gaus- 
sian mode. Therefore, the orbital angular momentum is 
not conserved in the parametric down conversion process, 
and the crystal is expected to suffer a twisting torque. 
This effect is analogous to the mechanical torque applied 
to a cjuater waveplate used for light polarization conver- 
fcyl: 

verse mode conversion |2j. 



sion [1|, or to a pair of cylindrical lenses used for trans- 



Thc reason for this assimetry in the OAM conservation 
of the pump can be explained when the propagation of 
paraxial beams in anisotropic media is investigated and 
the effect of cavity mode selection is considered. 



where ko — lu/c is the wave number in vacuum corre- 
sponding to frequency u>, and a is a constant to be con- 
veniently chosen. This constant will significantly sim- 
plify the paraxial propagation analysis in the birrefrin- 
gent medium. The constitutive relation D = e • E de- 
pends on the dieletric tensor e, that is diagonal when we 
use the crystalographic coordinates 



e — 
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ni 



(2) 



The wave equations for the electric field components 
can be derived from Eq.Q by using the constitutive re- 
lation and choosing a — to obtain 



4 d 2 x E x + d 2 y E x + d 2 E x - 1 - ^| d x d z E. 



kl n 2 x E x = 0, 



d x E y 



d l E v 



dtEy ~ ( 1 - -f ) dyd X E X 



) d v d z E z + k 2 n 2 E y 



(3) 



2 X E Z 



d 2 E z 



• d 2 E z 



d z d x E x 



k 2 n 2 z E z = Q 



Notice that for an uniaxial crystal (n x = n y ) we recover 
the equations obtained in ref. 29]. Let us now consider 
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propagation along an axis x' lying on the xy plane with 
an angle <f> with respect to the crystalographic axis x, 
as shown in Fig0 This definition of has the advan- 




z = z' 



which is the projection of the well known index ellipsoid 
on the xy plane. By making k x — nk cos <fi and k y — 
nko sin</>, we get 

cos 2 4> sin 2 </> 1 

2 1 T~ = ~2 W 

n x 

where n is the index of refraction for propagation along 
x' . On the other hand, if we substitute Eq.© in the 
constitutive relation D = e ■ E and in V ■ D = we find 
that 

n\ k x E x + n 2 y k y E y = . (9) 

Since n x ^ n y , this means that E and k are not orthogo- 
nal. Therefore, the Poynting vector S, that is orthogonal 
to E, is not parallel to k. Let us call <\>' the angle be- 
tween S and k . A straightforward geometric analysis 
allows one to obtain a simple relation between <fi and <fi': 



FIG. 5: Coordinate system used to describe the paraxial prop- 
agation through the anisotropic crystal. The walk-off angle 
(/>' is indicated as the angle between the Poynting vector S 
and the propagation axis x 1 . 

tage to match the angle and axis definitions usually given 
by crystal manufacturers. For our KTP crystal, cut for 
type II phase matching of 532 and 1064 nm, we have 
<t> = 23.5°. A rotated reference frame (x',y',z) can be 
used to describe the propagation inside the crystal. The 
coordinate transformation between the two frames is 



tan 0' = 



sin0 cos0(n 2 — n x ) 



(10) 



This angle is represented in Fig[S| It is related to the well 
known walk-off effect, which appears as a consequence of 
the crystal anisotropy. However, as we shall see shortly, 
the z polarized field will also have an anisotropic effect 
when the propagation of a transversely finite beam is con- 
sidered. This effect appears as an astigmatic deformation 
of the beam during the propagation along the crystal. 



x — x cos <p 
y = —x sin < 



y sin <p 
+ y cos < 



(4) 



z 



z . 



A. Plane wave analysis 

Two orthogonally polarized plane wave solutions prop- 
agating along x' can be found for Eqs.©. One with 
E x = E y = and E z ^ (z polarized) and another 
polarized in the xy plane {E z = 0). For the z polarized 
solution, only the last of Eqs.@ remain and its solution 
is 



E z = E 0z e 



in z ko x 



(•5) 



The plane wave solution polarized on the xy plane can 
be found by making E z = in the first two of Eqs.©. 
This solution is of the kind 



E = E e* ( fcx x+ky ^ , 



(6) 



where Eo = Eq x x + Eo y y. From substitution of Eq.JBJ 
in the first one of Eqs.jSl), we find 



2 p 

e i _y_ — i,2 

2 + 2 — Kl 



V 



(7) 



B. Paraxial propagation 

On ther other hand, to obtain a direct solution of eq.Q 
for a paraxial beam propagation is not so straightforward 
and some carefull aproximations have to be made to un- 
couple the diferential equations for each polarization. For 
the z component, the wave equation has the form 



E z - 



a - 



- d y d z E y 



. 



d x d z E x 



(11) 



To reduce this equation to the paraxial wave equation 
for z polarization, we can begin by eliminating the terms 
with cross derivatives. One way to do this is to approxi- 
mate the biaxial crystal by an uniaxial one for the z po- 



larization. This is valid since 



chose a 



we have 



>y\ 



= 10~ 2 for i = 
In the limit n. 



X 



-n y \. If we 

{x,y}, giv- 
ing a very small contribution, 
these terms will vanish, and we have the uniaxial crystal 
studied in ref . [29^ . 

A paraxial solution E z = u z (x' ,y' , z)e~ m * k " x of 
Ea. (|ll|) can be obtained if we adopt the rotated refer- 
ence frame. The resulting equation is close to the parax- 
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ial wave equation, except for the assimetry in the coeffi- 
cients of the transverse second-order derivatives: 



-4 d z u z = 2in z k d x >u z 



(12) 



The assimetry between the transverse coordinates y' 
and z appears as a rescaling of the z coordinate. This 
means that the optical beam follows an astigmatic prop- 
agation inside the crystal with diferent diffraction scales 
for each transverse coordinate. Let us separate the 
dependence of u z on y' and z making u z (x' ,y' , z) = 
U z {x' , y') V z {x' , z), in order to obtain two paraxial wave 
equations for the beam diffraction in each transverse di- 
rection: 



2in z k d x <U z 



-^d 2 z V z = 2in z k d w ,V z 



(13) 



When calculating the propagation of the beam through 
an OPO cavity, this diffraction assimetry can be seen as a 
different effective length of the crystal for each transverse 
dependence of the mode function. For a crystal with 
length £, the effective length for the U z propagation will 
be £/n z , as usual in the treatment of beam propagation 
through an uniform crystal |2^,|3^|. For V z , the effec- 
tive length it will be £n z /n 2 , resulting in an assimetry in 
the effective cavity length for each transverse evolution. 
The calculation of the cavity geometry, and the resulting 
beam parameters expressed by the Rayleigh length Xr, 
will therefore differ for the two transverse coordinates. 

Let us now turn to the paraxial solution for the field 
polarized on the xy plane. Since a plane wave solution 
with E z = can be found, it is natural to concieve a 
paraxial solution for which E z is negligeable. Therefore, 
if we choose a = n 2 in Eq. (JJJ and use the rotated coordi- 
nates, we obtain the following propagation equation for 
E y 

[cos 2 4> d 2 , + sin 2 4> d 2 , — sin 2<j> d x >d y i 

2 

Ti 

+-f (sin 2 0fi£, + cos 2 <t>d 2 y , + wa2</>d x .d y .) 



d 2 + k 2 n y ]E y =0 



(14) 



We now try a paraxial solution of the form E y = 
u y (x' , y' , z)e~ lnk ° x in Ea. lfT^ . using Eq.JSJ) and making 
the paraxial approximation to obtain: 



1/2 



2ih 



Qfly 



cos 



V ■ 2 

-4 sin - 

T 2 



[d X 'Uy + tail^'dy'Uy) = 



-| cos 2 , 



dy.Uy 



d 2 u v 



(15) 



where 0' is the walk-off angle given by Ea. (|10|) . In or- 
der to obtain a paraxial wave equation, a second coor- 
dinate transformation y" = y' — tan^'rr/ is necessary. 



This transformation corresponds to a transverse offset 
of the xy polarized beam. Using Eq.JSJ and defining 
£ 2 = sin 2 (j) + ( n l/n x ) cos 2 4> , we can rewrite Ea. lflTjl as 



2ik — d x >u y = £ 2 dy,,u y 



d 2 u 



v > 



(16) 



that is, the usual paraxial equation with wave vector 
ko(n 2 /n) and a rescaled transverse coordinate y"/£. 



However, since n x 



this transverse rescalling is 



much smaller than the one present in the z polarized 
field. Therefore, while the z polarization has a significant 
astigmatism but no walk-off, the xy polarization presents 
walk-off and a small astigmatism. From now on we shall 
designate the xy polarized field as the extraordinary wave 
and the z polarized field as the ordinary wave. 

As we made for the z component, we can try a factor- 
ized solution of the paraxial wave equation (| 1 fc>|) of the 
form u y — U v (x', y") V y (x' , z) , so that 

ed 2 „U y = 2ik (n 2 y /n)d x ,U y , 



d 2 V y = 2ik (n 2 y /n)d x ,Vy 



(17) 



Thus, the paraxial propagation inside the crystal is well 
described by Eas. l|13|) and i|17|) for the ordinary and ex- 
traordinary waves respectively. A paraxial equation for 
the x component of the extraordinary wave can be ob- 
tained on the same lines leading to Eq. l|17fl . 



IV. ASTIGMATIC CAVITY 

All scalling parameters appearing in Eas. (|13l) and H17|) 
can be absorbed by a suitable definition of an effective 
wave number fc e // for each transverse direction and for 
each polarization. This brings Eas. (|13fl and (|17fl to the 
general form 



d 2 yU{x,y) = 2ik eff d x U{x,y) 
The normalized solution of Eq. (|18|) is 



(18) 



U(x,y) 



k e ff XR 



1/4 



keff XR 

y\ i — — i x 



exp 



2(x + ixr) 



— I a ret an 



(19) 



where xr is the Rayleigh length, and H n (x) is the Her- 
mite polinomial of order n > 0. The term arctana;/xij is 
the well known Gouy phase shift. This term avoids mul- 
tiple resonances of high order Hcrmitc-Gaussian (HG) 
modes in a high finesse cavity for the signal and idler 
modes of the OPO. The beam propagation is character- 
ized by the beam waist wq — y / 2xR/k e ff and the wave- 
front curvature R(x) — x (l + x R /x 2 ) . The change in 
the effective wave number is equivalent (in terms of beam 
diffraction) to the propagation in a shorter length of free 
space. Since the effective wave number depends both on 
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polarization and transverse direction, we can consider a 
different propagation length in each case. 

Let us now consider the refractive index of the KTP 
crystal at 1064 nm (n x = 1.7404, n y = 1.7479, n z = 
1.8296) and 532 nm (n x = 1.7797, n y = 1.7897, n z = 
1.8877), according to the manufacturer, Cristal Laser 
S.A. We have, for the extraordinary wave, a refractive in- 
dex n(1064nm) = 1.7467 and n(532nm) = 1.7881. From 
the distance Lq = 17.4mm between the mirrors in our 
near-confocal cavity, and the crystal length £ = 10.0mm, 
we can calculate the effective length L of the cavity for 
each transverse mode, and for each polarization, in the 
infrared case. Using the relation 

L = L -£ kef ? - k ° (20) 

k eff 

we obtain, from the values of k e tf given by Eas ll3ll7l 

L° y , = 12.87mm, L° z = 13.40mm, 

Ly„ = 13.17mm, L e z — 13.12mm , (21) 

where the superscript o(e) refer to the ordinary (extraor- 
dinary) wave. The effect of the walk off for the extraordi- 
nary wave has been taken into account, but the correction 
was ~ 10" 4 , and could be neglected. The values of the 
Rayleigh length inside the cavity (x 2 R = L 2 (2R m —L)/4) 
for each transverse direction of the beam, and for each 
polarization, differ by less than 1%, and cannot be no- 
ticed in the free-propagating beam 

x° Ry i — 6.500mm, x° Rz — 6.497mm, 

x Ry „ = 6.499mm, x e Rz = 6.500mm . (22) 

On the other hand, the total Gouy phase shift ac- 
cumulated in a round trip inside the cavity $ = 
4arctan(-\/L/-v / 2i? m — L) will be 

= 3.122rad, $° = 3.204rad, 
<S>1„ = 3.167rad, §\ = 3.161rad. (23) 

The phase added in a round trip depends on the order 
of the Hermite-Gauss TEM mn mode resonating inside 
the cavity. The total Gouy phase for this mode is 

$ = (m + 1/2)$, + (n + l/2)$„. (24) 

From the calculated values of the Gouy phase shift, we 
see that there will be a small phase diference between 
the TEM i and the TEMi modes. This difference will 
result in a splitting of the resonance position. At 1064 
nm, this separation is of 82mrad for the ordinary wave, 
and 6mrad for the extraordinary one. 

In order to study this splitting, we pumped the OPO 
with a Laguerre-Gaussian (LG) mode obtained with an 
astigmatic mode converter 0. The LG mode is the su- 
perposition of two HG modes orthogonally oriented, that 
is, a TEMoi and a TEMi mode. Once the OPO cav- 
ity is scanned, a single resonance peak is expected if the 



cavity is degenerate for the two TEM modes. Otherwise, 
two resonance peaks are expected, one corresponding to 
each TEM mode. In FigHj] we show the resonance peak 
of a high finesse cavity for a 532 nm LG pump. The 
polarization of the pump laser was rotated in order to 
provide both, the ordinary and extraordinary waves. For 
the vertical polarization (ordinary wave), a double res- 
onance is observed as expected. This splitting shows 
a round trip phase difference of 88mrad, in reasonable 
good agreement with the predicted 93mrad for 532 nm. 
On the other hand, for horizontal polarization (extraordi- 
nary wave) , the LG resonance presents a single peak. In 
this case, the splitting is expected to be around 8mrad, 
well below the resolution of the cavity used for this mea- 
surement. 
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FIG. 6: Cavity resonance peak for a pump LG beam with 
a) horizontal (extraordinary) polarization and b) vertical (or- 
dinary) polarization. In the second case the resonance peak 
splits in two, clearly showing the simetry breaking between 
the two HG components of the LG beam. 

From this analysis we conclude that the OPO can sup- 
port the oscillation of an LG mode for the extraordi- 
nary wave, since its HG components have a degenerate 
(or quasi-degenerate) resonance frequency. On the other 
hand, an LG mode in the ordinary wave cannot oper- 
ate because its HG components will not have the same 
resonance frequency. This explains the results shown in 
Fig^J that is, the orbital angular momentum (OAM) is 
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transferred from the pump laser (extraordinary wave) to 
the idler mode (extraordinary wave) but not to the sig- 
nal mode (ordinary wave). Notice that, under our ex- 
perimental conditions, only one of the down converted 
modes oscillates in a high order transverse mode, while 
the other one oscillates in the fundamental transverse 
mode. So, the OAM exchange between pump, signal and 
idler modes is governed by the cavity dynamics under the 
crystal anisotropy, involving polarization and transverse 
profile aspects. 

V. THEORETICAL MODEL 

Transverse multimode operation of OPOs has already 
been theoretically discussed in Ref . |33j | . The pump beam 
can excite many different cavity modes for signal and 
idler, but in general it is the one with the lowest thresh- 
old that survives. Therefore, modes with best recover- 
ing integral should oscillate. To extend this description 
to our experiment, we must take into account the walk- 
off and the astigmatism due to the crystal anisotropy. 
As we have seen, the astigmatism will introduce a phase 
shift between the two Hermite-Gauss components of the 
Laguerre-Gauss beam. We can choose to treat the prob- 
lem either in the Laguerre-Gauss basis or in the Hermite- 
Gauss one. For the Laguerre-Gauss basis, the astigma- 
tism couples the right-handed beam to the left-handed 
one. In the Hermite-Gauss basis, this coupling implies 
in a phase difference between the two first order modes. 
Here we chose to work in the Hermite-Gauss basis, but 
the change of basis is straightforward. 

In order to study the dynamics of the relevant trans- 
verse modes, we shall consider the normalized mode func- 
tions Ujk{x' , y', z), where j — p,s,i for pump, signal and 
idler respectively, and fc = 0, h, v for the Hermite-Gauss 
TEM 00 , TEMio and TEM i respectively. The overlap 
integrals, 

A k i m = J J J u pk u* sl u* m dx' dy' dz , (25) 

play an important role in the dynamics since they deter- 
mine the transverse mode coupling. The mode functions 
Ujk(x',y',z) are given by Ecis. lfH^ . l(T7|l and (|T§1) . where 
astigmatism and walk-off effects are taken into account. 
The walk-off is slightly different for pump (4.1 mrad) and 
idler (3.2 mrad), and the significant astigmatism occurs 
in the z direction of the signal mode. The integrals are 
calculated in the whole crystal volume. 

With the overlap integrals, we can obtain the dynamic 
equations for the transverse mode amplitudes. From all 
possible combinations of oscillating modes, the cavity 
parity will restrict the number of transverse modes for 
a given longitudional mode. If there were no anisotropic 
effects, with a first order Laguerre-Gaussian pump mode, 
which is odd, signal and idler must have oposite parities 
in order to give a nonvanishing overlap integral. There- 
fore, for isotropic propagation, if signal oscillates in a 



first order mode, idler must oscillate in the fundamen- 
tal one, and vice-versa. In principle, this parity selection 
brakes down for an anisotropic medium specially due to 
walk-off. However, when the overlap integrals are calcu- 
lated, we can see that the integrals for odd combinations 
of modes, like (v, 0, 0) or (v, v, v) for example, are indeed 
much smaller than those obtained with an even combina- 
tion like (v, v, 0). This allows us to neglect many of the 
mode couplings and restrict the number of dynamic equa- 
tions. Two kinds of operation regimes are observed: ei- 
ther the signal beam oscillates in the fundamental TEMoo 
mode, while the idler lies in the TEMoi and io subspace 
(peaks 1, 2 and 4 in Fig[3J), or the idler beam oscillates 
in the TEMoo mode (peak 3 in Fig|3|. Let us describe 
these regimes separately. 

A. Signal beam operating in the TEMoo mode 

In this case, the set of dynamic equations for pump, 
signal and idler transverse mode amplitudes is: 

a pv = -[t p + «(A P + <T p )]a pv - ixKov a soa lv + E in /V2 

a p h = ~h P + «(A P - cr p )]a ph - ixA* hQh a s0 a ih - i E ir J a/2 

a s o = -(7 + *A s )a s0 + ixKo v a pv a* v + ixkhoha ph a* h 

a lv = -[7 + i(Aj + (Ti)]a iv + ixKova pv a* s0 

dih = -[7 + - ai)]aih + iX-kh.oh.aph.alo- ( 26 ) 

where the subindices p, s and i refer to pump, signal 
and idler respectively, and 0, v and h refer to fundamen- 
tal (TEM o), vertical (TEM i) and horizontal (TEMi ) 
transverse modes. Pump losses are described by 7 p while 
a common decay rate 7 represents the losses for signal 
and idler. The respective cavity detunings for pump, 
signal and idler are A p , A s and Aj. The astigmatic sym- 
metry breaking is accounted for through the frequency 
splitting parameters a p for pump and <Ji for idler. They 
are calculated with the help of Eq. I|24|) . The pump beam 
amplitude transmitted through the input mirror is rep- 
resented by the source term Ei n . Since it is prepared in 
a Laguerre-Gauss mode, the source terms appearing in 
the dynamic equations for the amplitudes a pv and a p h 
are tt/2 out of phase. Finally, x is the nonlinear coupling 
constant. 

The dynamic equations, as well as their steady state 
solutions, are considerably simplified if we express time 
in units of the cavity round trip time r and define the 
following normalized variables 

bjk = X Aooo r ajk , x in =X Aooo r E in , 
7j = 7i 7- , Aj = Aj t , <jj =Ajt, (27) 
_ Afe 1 m 

T]k Ira "7 • 

A000 

As before, j = p,s,i for pump, signal and idler respec- 
tively, and each of the subindices k, I and m may assume 
the values 0, h or v. Cavity losses are around 29% at 
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532nm and 1% at 1064nm which gives j p = lA5mrad 
and 7 = 5 mrad . In the absence of astigmatism and 
walk-off the relevant normalized overlap integrals are 
Vwo = Vhho = r/vOv = rjhoh = 0.71 approximately. 
When the walk-off effect is considered, the overlap in- 
tegrals are averaged over the crystal volume. More- 
over, the astigmatism is included through the appropri- 
ate correction of the mode functions. Taking into ac- 
count the experimental values for the walk-off and astig- 
matism parameters we find ri vv0 — 0.70, %/i0 = 0.60 
and r) v o v ~ %o/t ~ 0.71. So, a significant change is 
obtained only for r/h ho- 
lt is instructive to consider the steady state solution 
of Eas. (|26|) in the simplified condition A s = A» = a p — 
<ji = and rjvov — r/hoh — V which correspond to neglect- 
ing walk-off and astigmatism. In this case the orbital 
angular momentum is perfectly transferred to the idler 
beam which will also oscillate in a Laguerre-Gauss mode 
with the same topological charge of the pump beam. 
Therefore, the steady state solutions are 



h- = , I p 



I s q — h-\ 



h 



7 2 



7 P 



,(28) 



-2- 

tn 
c 

B 




t(T) 



FIG. 7: Time evolution of the down converted beam inten- 
sities (in units of the cavity round trip time) obtained from 
numerical integration of the dynamic equations 1261 . The in- 
set shows the corresponding expected images for signal (s) 
and idler (i). The parameter values used are A p = 0.071 y p , 
A s — , Ai — 1 mrad , 7 P = 145 mrad ,7 = 5 mrad , a v = 
4 mrad , Oi = 3 mrad , Xi„ = 3 il and i),o« = r\u o h = 0.71 . 
The horizontal solid line shows the analytical value To- 



where we defined the normalized intensities Ijk — \bjk\ ■ 
The Laguerre-Gauss amplitudes bj± are given in terms 
of the Hcrmite-Gauss amplitudes as 



b jv ± i bjh 
V2 



(29) 



The threshold value of Xi n for parametric oscillation is 
obtained by setting Iq = so that 



XL 



rf 



(% 2 + A^) 



(30) 



As we shall see, a different threshold condition is obtained 
for the other operation regime, in which the idler beam 
operates in the TEMoo mode. 

The analytical solution for the steady state includ- 
ing all parameters is cumbersome but Eas. (|28|l give us 
a good estimate for the orders of magnitude. In fact, 
as we discussed in section IV, the expected value for the 
pump and idler splitting parameters are ideed very small, 
dp = 4 mrad and &i — 3 mrad (the splitting parameter is 
half the astigmatic phase shift calculated from Ea. (|24|l h 
However, this small splitting may be responsible for par- 
tial transfer of the orbital angular momentum from the 
pump to the idler mode. In order to illustrate this, we 
numerically integrated the dynamic equations (|26|l with a 
fourth order Runge-Kutta method until the steady state 
was reached. In Fig[7| this time evolution is shown to- 
gether with the value Io given by Ea. l(2%|) . In the inset, 
we show the expected image for signal (S) and idler (I) 
obtained with the numerical steady state results. A good 
qualitative agreement is obtained with the experimental 
results corresponding to peaks 1, 2 and 4 of FigEI 



B. Idler beam operating in the TEMoo mode 



In this case the dynamic equations are: 



Qpv — 


-hp 


Fi(A p 


+ &p)]a P v 


- ixKvO a sv a M + E in/V2 


a ph — 


-hp 


F i(A p 


- Vp)]a P h 


- ixK,ha a shaio - i E in / \[2 


(Lsv — 


-[ 7 4 


-i(A s - 


1- <r s )]a sv 4 


iX^woa pv a* (31) 


d s h = 


-[74 


-i(A s - 


- <r s )]a sh 4 


- iX^-hhoa p ha* 


dio — 


-(7 4 


- iAi)a 


m + ix^-vvO a pva* sv 4 ixh-hhoa-phO* s hi 



The transverse mode splitting now appears in the dy- 
namic equation for the signal beam and is represented 
by the parameter a s . However, the splitting parameter 
is expected to be of the order of 41 mrad. Since cavity 
losses in the infrared are of the order of 1%, the corre- 
sponding normalized decay rate is 7 = 5 mrad, so that 
<7 S >> 7. Under such conditions it is impossible for 
the OPO to support the simultaneous operation of the 
h and v modes necessary to compose a Laguerre-Gauss 
mode. Therefore, the orbital angular momentum cannot 
be transferred to the down-converted beams. The cavity 
tuning will select the signal Hermite-Gauss mode whose 
resonance frequency is closer to the idler resonance. For 
example, for A s = —o s the cavity frequency falls far 
away from the h signal resonance while the v mode gets 
on resonance. In this case a s h ~ and the steady state 
solution of Eas. (|31J) can be analytically obtained. Notice 
that the normalized overlap integral rjhho will not play 
any role in this case. We therefore set r) = T) vv q and use 
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the same normalizations adopted in Eas. (|27(l to find 



0. 



IiO — Isv — — ~~" 




7 2 



(32) 



The h component of the pump beam does not couple 
do the down-converted modes so that its steady state 
solution corresponds just to an empty cavity. On the 
other hand, the v component of the signal beam, as well 
as the fundamental idler mode, presents a steady state 
intensity lower than the one found in Eas. H28J) for the 
same pump level Xi n . This corresponds to the situation 
found in peak 3 of Fig|3 which is clearly lower than the 
other infrared peaks. Again, the oscillation threshold is 
readily obtained by taking I' G = 0: 



XL 



^! (f + A 2 



(33) 



It is twice the threshold value for the case where the or- 
bital angular momentum is transferred for the idler beam, 
what is also coherent with the lower height of peak 3 in 
FigEl 

The numerical evolution of Eas. (|31l) using a fourth 
order Runge-Kutta method was performed without the 
simplfying assumptions. These results are presented in 
FiglHl where the inset shows the expected images for sig- 
nal and idler. The walk off and astigmatic effects were 
fully considered and a good agreement with the exper- 
imental result was obtained. The theoretical model de- 
veloped here will be useful for future investigation of the 
transverse mode dynamics in the quantum domain. In- 
teresting perspectives can be envisaged if the OPO oper- 
ation is subject to an injected signal. Recent studies on 
degenerate |3jj and nondegenerate parametric pro- 
cesses with injected signal has considered interesting is- 
sues like the preparation of quantum correlated states 
(Einstein-Podolsky- Rosen states) as well as the study of 
critical behaviours of the OPO operation |3S|. 



1 

0.9 

0.8 

0.7 

J_0.6 

io.5 
tn 

I 0.4 
c 

0.3 
0.2 
0.1 






2000 3000 

t(x) 



FIG. 8: Time evolution of the down converted beam inten- 
sities (in units of the cavity round trip time) obtained from 
numerical integration of the dynamic equations 1311 . The in- 
set shows the corresponding expected images for signal (s) 
and idler (i). The parameter values used are A p = 0.28 7 P , 
A s = —41 mrad , Aj = , 7 P = 145 mrad ,7 = 5 mrad , 
<7 P = 4 mrad , a s — 41 mrad , Xi n = 1.5 xl and rj v v = 0.70 . 
The horizontal solid line shows the analytical value Iq. 



Our results open interesting perspectives related to the 
quantum regime of the OPO operation. The experimen- 
tal setup can be used to generate twin beams (intensity 
difference squeezing) and the production of the optical 
vortices under squeezed operation can lead to interesting 
quantum patterns in the down-converted beams. More- 
over, the quantum dynamics of the type II operation 
may lead to entanglement between the polarization and 
the spatial degrees of freedom [3jJ . We are theoretically 
investigating this possibility and an experiment will be 
setup in the near future. 
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VI. CONCLUSION 

We have shown that the transfer of orbital angular 
momentum in intracavity parametric down-conversion is 
strongly subjected to cavity and anisotropy effects. This 
can be achieved if pump, signal and idler are in a set 
of modes where the Hermite-Gauss components of the 
Laguerre-Gauss mode are degenerate inside the cavity. 
While that can be easily achieved for the idler beam, the 
signal beam cannot fullfill this condition unless cavity 
losses are large. The signal beam can still oscillate in a 
higher order transverse mode, but in this case, the thresh- 
old power increases, and the orbital angular momentum 
is not transferred to the outcoming beams. 
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